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Введение. При численном решении задач теории упругости в трехмерной постановке методом конечных 
элементов применяются конечные элементы (КЭ) в форме параллелепипедов, призм и тетраэдров. Обычно 
построение матриц жесткости объемных КЭ базируется на принципе изопараметричности, суть которого 
состоит в использовании для аппроксимации геометрии и перемещений полиномов Лагранжа. В расчетной 
практике наибольшее распространение получили так называемые полилинейные изопараметрические КЭ 
с линейным законом аппроксимации перемещений. Главный недостаток данных элементов кроется в эффекте 
«юос_те» («запирания») при моделировании изгибных деформаций. Причем погрешность численного решения 
существенно возрастает в случае, когда конструкция, по сравнению с обычными деформациями, претерпевает 
значительные смещения как жесткое целое. Многолетний опыт решения задач механики деформируемого 
твердого тела методом конечных элементов показал, что существующие объемные КЭ обладают медленной 
сходимостью при моделировании изгибных деформаций пластин и оболочек. Цель настоящего исследования 
состоит в построении на основе метода двойной аппроксимации матриц жесткости полилинейных объемных 
КЭ повышенной точности, позволяющих учитывать жесткие смещения. 

Материалы и методы. Для построения матриц жесткости объемных КЭ применен математический аппарат 
метода двойной аппроксимации, суть которого состоит в раздельном представлении функций распределения 
перемещений и деформаций внутри элемента. Хранение и обработка результирующей системы уравнений 
реализованы в алгоритмических терминах разреженных матриц. Разработка программного обеспечения 
и проведение вычислительных экспериментов осуществлены с использованием 64-х разрядной вычислительной 
платформы Мисгозой У15иа1 Зо 2013 и компилятора пе® РагаПе!| Зию ХЕ 2019 со встроенным текстовым 
редактором Пие!® У15иа1 Еогпгап Сотрозег ХЕ 2019. Визуализация результатов расчетов выполнена с помощью 
дескрипторной графики пакета компьютерной математики МаЙаб. В качестве тестового образца использован 
объемный восьмиузловой КЭ $ОГЛО185 программного комплекса АМЗУ$ Месратса|. 

Результаты исследования. Разработано математическое и программное обеспечение для исследования 
напряженно-деформированного состояния массивных конструкций при различных видах внешнего 
воздействия. На тестовых примерах с известными аналитическими решениями выполнена верификация 
авторизированного пакета прикладных программ. Показано, что построенные КЭ по точности удовлетворяют 
основным требованиям, предъявляемым к конечно-элементному моделированию пространственных задач 
теории упругости. 

Обсуждение и заключение. Проведенное тестирование разработанного математического и программного 
обеспечения показало, что построенные на основе метода двойной аппроксимации конечные элементы успешно 
конкурируют с аналогичными объемными элементами $ОГ/ЛО185 программного комплекса АМЗУ$ Месвашса|. 
Предлагаемые элементы могут быть интегрированы в отечественные импортозамещающие программные 
комплексы, реализующие метод конечных элементов в форме метода перемещений. 


© Гайджуров П.П., Савельева Н.А., 2023 


Механика 


365 


Бйр://уезиК-доп$а.ги 


366 


Афуапсеа Епятеетия Кезеагсй (Коз1оу-оп-Ооп). 2023;23(4):365-375. е155М 2687—1653 


Ключевые слова: метод конечных элементов, моментная схема метода конечных элементов, метод двойной 
аппроксимации, объемные конечные элементы, тестирование конечных элементов 


Благодарности. Авторы выражают благодарность рецензентам, чья критическая оценка представленных 
материалов и высказанные предложения по их совершенствованию способствовали значительному повышению 
качества настоящей статьи. 


Для цитирования. Гайджуров П.П., Савельева Н.А. Применение метода двойной аппроксимации для построения 
матриц жесткости объемных конечных элементов. —Айуапсей Епетеетия Кезеагсй (Козюу-оп-Ооп). 
2023;23(4):365-375. № рз://401.0ге/10.23947/2687-1653-2023-23-4-365-375 


Опета[ агие 
АррИсабоп оЁ Фе РоипЫе Арргохипайоп Ме@од ог Соп$гисйпо ЗЯЙпе$$ Майлсе$ 
о? Уоштейле Ешце Еетепт$ 


Реег Р. Са147Вигоу ©. Мша А. бауееуа@ 
Доп Зае Теснтса1 Отуегзиу, Возюу-оп-Ооп, Кизз1ап Еедеганоп 


эрр-161@уапа4ех.га 


АЪзтгасе 

ттодиспоп. УПеп питепсаПу зоутя ргоепаз оЁ е!азйсИу ШФеогу ш а гее-4итепз1опа| Коглианоп Бу Фе йпце 
е!етепё пеоа, Нпие @етеп6 (ЕЕ) ш Фе Югп о# рагаПе!ерре4$, ризт$ ап4 {ейтаБейга аге изе4. Кезшайу, пе 
сопзгисНоп оЁ $#пез$ тайлсез оЁ уопитейлс ЕЁ 1$ Базе оп Фе рипср[е оЁ 1зорагатефеИу, усЬ шуоуез ве 
Гастапее ро[упоп1а1$ 0 арргохипайе фе эеотейу апа @1р1асетеп. ш сотрщайопа] ргасйсе, Фе п105{ \4езргеаа 
ЕЕ аге Фе зо-саПе4 ши Йтеаг 1зорагатейлсе ЕЁ \/ИВ а Ппеаг |а\у о{ арргохипайоп оё @1зр1асетеп. Тве тат 
Ч1за4уатиасе оЁ Шезе е1етепз Пез ш 1е ‘Чоскте” еЙес{ \Пеп подайте Беп4тз деЮгтаНнопз. Могеоуег, фе еггог 
ог Фе патепса| зоайоп шсгеазез Чгазиса!у ш Фе сазе \Веп Фе касте, ш сотрапзоп  сопуепйопа| 
еогтайопз, ип4егоое$ 51ет!Исап 415р]асетеп{5 аз а пол уБое. Гоп2-епп ехрепепсе ш 50|уше ргоепаз оЁ 
ЧеРогта Ме зоП4 тесвап1лс$ Бу Фе ВпИе е!етепё тефо@Я Баз зБо\’п Фа ех15бпе уотейле ЕЁ Бауе $10\ 
сопуегеепсе, зресИсаПу, увеп то4деПпе Беп4те АеюогтаНоп$ оЁ рез апа вез. ТЬ$ зеа4у ай аё сопзгасипе 
И ипез$ пайлсез оЁ пиИтеаг уо[итее ЕЕ оЁ шсгеазе4 ассигасу аПо\уше Юг 111А 91р]асетеп Базе оп е 


4очЫе арргохипаНоп тефо4. 

Миенаб апа Мешо45. ТВе тафетайса| аррагайаз оЁ фе дочЫе арргохипаНоп те#фо4 Ъазе4 оп Фе рипсф!е оГа 
зерагайе гергезепаНоп оЁ йе Ч151БаНоп Впсйопз оЁ а1зр!асетеп!$ ап деРогтаНоп$ шз14е Фе е]етепть \уаз изе4 о 
соп$гисЕ Фе зи пез$ тай1сез оР уоатейлс ЕЕ. Те зюогазе ап4 ргосеззте оЁ Фе тез пе зузет оЁ едиаНоп$ \’аз 
пир!етегие4 ш а1еогИс {1егпл$ оЁ зрагзе тайлсез. Зой\аге аеу@юортеп( ап@ сотрщайопа| ехрегилеш \уеге 
саглед оц изте фе Мтсгозой У\1за1 За 2013 64-65 сотрийпе р!аогт апа Фе Пе] ® РагаПе! Зо ХЕ 2019 
сотрИег УИ Фе ищеотаеа шие ® Ува! Еойгап Сотрозег ХЕ 2019 {ехё е4йог. У1зиаПтаНоп оЁ йе сасийайНоп 
тезиНз \уаз регРогте изше фе 4езсирюг этарШс$ оР фе МАТГАВ сотшрщег тафетайс$ расКазе. А 1агое еазЪ- 
по4е $011О185 СЕ ор \е АМЗУ$ МесВатса| зоЙ\аге сотр[ех Уаз изе4 аз а {е5{ затр/е. 

Кезий$. Мафетайса| {001 ап зоЙ\аге \еге 4еуе]оре4 ю з@а4у Ше зге$$-5гали зе оЁ таззтуе згасигез ипаег 
уапоц$ (урез о еж{егпа| асйоп$. ТБе аиогте аррИсайоп зоЙ\аге расКазе \’аз уемйе@ оп {е5ё ехатар|ез ИВ 
Кпоууп апа[уйса| зо оп$. Ц Ваз Бееп зВо\уп фай пе сопзиласе4 ЕЁ ассигайе]у зайзКу Фе Баз1с гедитетепт Рог Йпие 
е!етепё тодеНп? оЁ зрана1 ргоеплз оЁ е!азйсйу Феогу. 

Оёсизяюоп апа Сопияон. Тре реогте4 {езНпе оРФе деуеоре4 тафетайса] ап4 ргостат юоКИ Баз зпо\уп фай йе Япие 
еетеп сопзиасе4 оп Фе Ба$15 оЁ фе дочЫе арргохипаноп тейфо@ сап зассеззШу сотрее у зипЙаг $0185 
уо[лтейлс е]етеп5 оЁ Че АМЗУЗ$ МесКашса| зой\уаге сотр/ех. ТВе ргорозе4 @еглеп сап Бе пщеотае шю Чотезйс 
пирой-заб$ийте зойхуаге зузёетлз фа ппретеп( ве Япие @етеп плево4 т Фе Юпп оРе а1р!асетеп тефод. 


Кеууог45: НпИе еетепё тео, потепё зсВете ог Нпйе еетепё тефо4а, доче арргохипайоп тефо4, 
уо[итейлс Ипие е]етлеп{5, АпИе е]етепё {еп 


Аскпо\едоетеп 5. ТВе аифогз арргеслайе Фе геуле\уег$, \Позе сиса| аззеззтепе оЁ йе забиие таепа1$ ап4 
зиссезноп$ Бере4 {0 з1ютЙсапИу ппргоуе Фе диаШу о 15 агисе. 


Ког сИайоп. Са1А7ригоу РР, Зауееуа МА. АррПсайоп оЁ Фе оч е Арргохипайоп Мефоа ог СопзгасИпе 
ЗИРпез$ Майлсез оЁ Уоштешмсес Ешие Еетеп. Айуапсеа Епетеетие Кезеагсй (Коз1оу-оп-Ооп). 
2023;23(4):365—375. 6 рз://401.0г2/10.23947/2687-1653-2023-23-4-365-375 
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Введение. При конечно-элементном моделировании напряженно-деформированного состояния массивных 
тел используются объемные конечные элементы (КЭ) в форме параллелепипедов (гексаэдров), призм 
и тетраэдров, построение матриц жесткости которых выполнено, как правило, по изопараметрической 
технологии [1-5]. Вместе с тем известно, что полилинейные изопараметрические КЭ при использовании 
однослойной схемы неудовлетворительно моделируют изгибные деформации даже при существенном 
сгущении сетки [6, 7]. Суть данной проблемы состоит в эффекте «/юс®ие» («запирания») элемента вследствие 
так называемой деформации «ложного сдвига» [8, 9]. Для «улучшения» изопараметрических КЭ используют 
аппарат несовместных элементов, созданных путем введения дополнительных неузловых степеней свободы или 
вспомогательных аппроксимирующих полиномов [8]. Вместе с тем наиболее эффективным способом решения 
проблемы «заклинивания» КЭ является применение моментной схемы метода конечных элементов, 
теоретические основы которой были разработаны А.С. Сахаровым [7]. В последующем данный подход получил 
название метод двойной аппроксимации (МДА) [6]. Концептуально МДА базируется на раздельном 
представлении функций распределения перемещений и деформаций внутри элемента. Целью настоящего 
исследования является построение на базе МДА и тестирование новых объемных полилинейных КЭ, 
позволяющих моделировать поведение различных конструкций при различных видах внешнего воздействия. 

Материалы и методы. Рассмотрим семейство объемных КЭ, состоящее из восьмиузлового 


и шестиузлового элементов, в глобальных декартовых осях 2,, т=1,2,3 (рис. 1). Геометрию 


т 


и перемещения КЭ представим в следующем виде: 


п. п", 
= (К) р = (К) 
2 =>. = ФСО, х,,х,) Ия =У и ФС х,,х,), 
К=1 


К=1 
ую 


т ›? т 


где 2 — узловые координаты и перемещения; ф,(х,,х.,х.) — «функции формы», представляющие 
собой произведение одномерных полиномов Лагранжа первой степени; х,, х›, х, — локальные в общем 
случае неортогональные координаты КЭ; и, — число узлов элемента. Для базового восьмиузлового элемента 


(рис. 1 а) п, = 8 «функции формы» определяем по формуле: 
1 3 

Фи хх) = Ч+рих,), (1) 
г=1 


здесь р,„, — координаты узлов в локальных осях. Значения Р,‚ задаем в форме матрицы: 


РК 


РЕ ЕЕ Е 


7 


ИА 
Рис. 1. Объемные КЭ: а — восьмиузловой; 6 — шестиузловой 


Зависимость между ковариантными компонентами тензора деформаций в локальном базисе и 
перемещениями в глобальных осях имеет вид [7]: 


1 
5; == ( 2т,] ии: 2: а 
2 
(суммирование по повторяющемуся индексу) 
где 2„,=02„/дх,; и„,=ди„/дх,, 5=Ь]. 
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Связь между вектором деформаций {=} и вектором узловых перемещений {и} представим в матричной 


форме: 
{}=ФИ из, 


где блочная матрица [2] = [2 [Р1,..[В,, | 


(6хи,) 
субматрица 
[В], =[{2°}, (ро, 12°] ‚К =12,...п,. 


(6х3) 
Выражения для вектор-столбцов рассматриваемых КЭ имеют следующий вид [10]: восьмиузловой элемент 
(рис. Та): 


Ри |2 ил + (2и12 Е Рик )хо + (2илз +21 Ре) Хз г 


+(2 „23 +212 Рзк 12 из РЖ +21 Рок Ри) Хх Хз | 


Рэк | т? о +22 Ри )х, на +2 ш,2 Рзк)хз + 


+2 ила +2 12 Рзк 12 из Ри 122 Ри Рзь)х! Х: | 


Рзк |2 из (в +2 из Ри )х а +2 из Рок )х, + 


1 Ц 
(т) р 25 о р 5 я 
ни #23 а Ри) хи Х. | ь (2) 


211 Рок 22 Ри + (2 из Рэь 1 Рок Рук из Ри + 


„Ри Ри) 


211 Рук т2из Ри нь Рак +21 Рок Рк + 2тз Ри + 


+2 3 Ри Рэк)х, 


2 то Рук Т2тз Рок ЕР +22 Рик Рак +2 из Ро ыы 


+2 из Рак Рэь)х 


шестиузловой элемент (рис. 1 6): 


| 
з 
5: 
12 
т 
+ 
— 
м 
15 
=) 
С 
+ 
| 
з 
5: 
м 
| 
— 
< 
| 


2 из Физ + (Рив Физ + Раз ф 3 ре 
1 о Фа Ра Фен ) 


{р М ЕЕ - 
, = а Е ыы 
8 | 21 Фл2 +22 Фу: + (из Фи ил Физ + 


Здесь введены обозначения: 


г _ 92 ое О’я. 
та > бтаВ — 5 
х 
9х. ХЕХ =хз=0 дх.д в ХЕХ =Хз= 
ы 9°2„ 
т,123 
дх, дх, дх 
1 й з жехо=хз=0 
- _9ф 9 а 
а ав 5 Ох ; о, В =Ъ2,3. 
о. © В ХЕХ =Хз=0 


ХЕХ =Х. =0 
Выражения для «функций формы» Ф,(х,х,,х,) шестиузлового КЭ, полученные на основании 


полинома (1) с помощью принципа «вырождения», имеют следующий вид: 


| 
а 


. 1 
Ф, о 
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й 1 
Ф: 5 [+ ры) + рых) + рззх,) + 1+ ра) + рых) + рых) 
ь 1 
Фф. Рае ыыы 
Е 1 
ф; Рю Ва) ; 


— 
ф. = [И+рух) + рых.) + рух) + + рых) + рых.) + рых.) 


Формулы (2) и (3) являются основой для построения матриц жесткости рассматриваемых КЭ. 
Соответствующее программное обеспечение разработано на базе вычислительной платформы Мисгозой У15па1 
Эшаю и компилятора П\е® РагаПе| Зо ХЕ с встроенным текстовым редактором Пие® У15йа| Еоггап 
Сотрозег ХЕ. Процессы хранения и обработки глобальной матрицы жесткости реализованы в терминах 
разреженных матриц [11]. Для визуализации результатов расчетов использована дескрипторная графика 
компьютерной системы МайаБ. 

Результаты исследования. Исследование точности и сходимости разработанного конечно-элементного 
алгоритма выполнили на тестовых примерах, имеющих аналитическое решение. В тестовых примерах 
приведены численные решения, полученные с использованием разработанных элементов и аналогичного по 
размерности элемента $0110185 программного комплекса АМЗУ$ Меспашса| [5, 11]. Ниже приведены 
примеры, подобранные таким образом, чтобы в них присутствовало сочетание изгибных деформаций и жестких 
смещений КЭ. 

Пример 1. Разрезное кольцо, жестко закрепленное в одном сечении и загруженное сосредоточенной силой 
на свободном конце. Расчетная схема кольца показана на рис. 2. Исходные данные: средний радиус КЁ = 0,2 м; 
размеры поперечного сечения 0,2х0,2 см; модуль упругости Е= 10" Н/м?; коэффициент Пуассона у = 0,3; 
сосредоточенная сила Р = 10 Н. 

Прогиб в точке приложения силы по теории кривых стержней (точное решение) составляет [7]: 


Ед В? 10.3,14.0,23 
= ЧЕ 


0,00302 м. 


ЕЛ 1.10'.8. 333-108 


Рис. 2. Расчетная схема разрезного кольца 


Результаты сходимости представлены в таблице 1. 
Таблипа 1 
Результаты сходимости для разрезного кольца 


1/6, м/% 
Сетка 
МДА АМЗУ$ 
2х2х32 —0,00170 / 44 —0,000433 / 86 
2х2х64 —0,00292 / 3,3 —0,00125 / 59 
2х2х128 —0,00293 /3 —0,00251 / 17 


В таблице 1 в числителе приведено значение прогиба в знаменателе — относительная погрешность б. 
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Визуализация перемещений и ‚ , полученная на базе МДА и АМ$У$ на сетке 2х2х64, представлена на рис. 3 


и 4. Отметим, что поле вертикальных перемещений, полученное с помощью АМЗУЗ, не отражает зоны 


с экстремальным значением и =-0,00314 м, показанной на рис. 3 стрелкой. 


и,,м 
-0.003138 
-0.002456 
^^ -0.001793 
-0.001121 
-0.0004484 


Рис. 3. Картина распределения И (МДА) 


и, м 
. 001346 
. 001186 
. 001026 
уе в х а . 865Е-03 
7 . 705Е-03 
. 545Е-03 
. 385-03 
.225Е-03 
. 646-04 
. 956-04 


Рис. 4. Картина распределения И, (АМЗУ$) 


® 
> 


В данном примере особенно отчетливо прослеживается важность учета жестких смещений. 

Пример 2. Квадратная пластина, жестко защемленная по контуру и нагруженная равномерно 
распределенной нагрузкой. Исходные данные: длина стороны а= 1 м; толщина й= 0,01 м; модуль упругости 
Е = 10$ Н/м?; коэффициент Пуассона у= 0,25. 

Точное значение прогиба в центре пластины определяется по формуле [12]: 

4 
94" 


и Б 


2 


ЕЮ 
тде а= 0,00126; Р= 12а) — цилиндрическая жесткость; 4= 0,00888889 Н/м? — интенсивность 
—у 
распределенной нагрузки. Точная величина и, (в метрах) равна коэффициенту %. 


В данном примере рассматривалась и часть пластины с учетом условий симметрии. Результаты сходимости 


в виде графиков И, —5 для однослойной и двухслойной моделей представлены на рис. 5 и 6. Здесь и далее под 


слоями подразумевается разбивка пластины на КЭ по толщине. На этих рисунках значения параметра 5 
1,2,3,4 соответствуют сеткам: 4х4, 8х8, 16х16, 32х32. Приведенные графики отображают результаты 
решения, полученные с помощью АМЗУЗ (линия 1) и МДА (линия 2). Горизонтальная линия, обозначенная 
цифрой 3, соответствует точному решению. 


Гайджуров П.П. и др. Применение метода двойной аппроксимации для построения матриц жесткости объемных конечных элементов 


Рис. 5. Графики И, -5; Рис. 6. Графики И, -5; 


однослойная схема двухслойная схема 


Из рис. 5 следует, что при однослойной схеме разбивки значения относительной погрешности на сетке 
32х32 составляют: $5ОГЛО185 — 16 %; МДА — 10,5 %. При использовании двухслойной схемы (рис. 6) на сетке 
32х32 имеем: 5011185 — 36 %; МДА — 2,8 %. 

Патч-тестирование КЭ выполнено для схемы разбивки 16х16х2 с искажением сетки (рис. 7). Результаты 


патч-теста в виде картин распределения прогибов и. приведены на рис. 8и9. 


Рис. 7. Схема разбивки пластины для патч-теста 
и ,м и ,м 
-0,001263 
су -. 958Е-04 
7 -. 851Е-04 -0,001010 
Е -. 745Е-04 
-. 639Е-04 -0,0007575 
И 
. 532Е-04 
СУ -. 426Е-04 -0,0005050 
| 
-. 319Е-04 
Е В ы -0,0002525 
. 213Е-04 
Г 
-. 106Е-04 
И 0,0 
Рис. 8. Распределение и, ($0Г10185) Рис. 9. Распределение и, (МДА) 


Как видно из рисунков, искажение сетки при использовании 30110185 приводит к более заметной 
асимметрии поля и, ‚ чем при использовании МДА. При этом значение максимального прогиба для КЭ МДА 


и, = 0,001263 м совпадает с точным решением. 


Пример 3. Круглая пластина, жестко защемленная по контуру и нагруженная равномерно распределенной 


нагрузкой. Радиус и толщина пластины: А = 1м; Й= 0,01 м. Механические константы аналогичны данным 


примера 2. 
Точное значение прогиба в центре пластины определяется по формуле [12]: 
_ 1 да‘ 
‘р 


При значении интенсивности распределенной нагрузки 4 = 0,00888889 Н/м? значение и,= 0,01563 м. 
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В процессе тестирования были использованы два варианта дискретизации 4 части пластины (сектора) на КЭ. 
В первом варианте три стороны сектора делились на равное число отрезков. Второй вариант базировался на 
радиальной регулярной схеме разбивки сектора. При этом число элементов вдоль радиуса и циркульной части сектора 
принималось одинаковым. Рассмотренные варианты дискретизации сектора для сетки 32х32 приведены на рис. 10. 


Вариант 1 Вариант 2 


хх 


Рис. 10. Варианты дискретизации сектора круглой пластины на КЭ: 1 — равномерная разбивка трех сторон сектора; 
2 — радиальная регулярная разбивка 


Результаты сходимости в виде графиков и, —5 для дискретизации 4 части пластины по варианту | при 
однослойной и двухслойной схемах разбивки представлены соответственно на рис. 11, 12 и рис. 13, 14. 


1 2 3 4 
|| И .9 1 2 3 45 


-0,002 -0,002 
-0,004 
-0,006 
-0,008 
-0,010 
-0,012 


-0,014 


-0,004 
-0,006 
-0,008 
-0,010 
-0,012 
-0,014 


и, М и, М 
Рис. 11. Графики и, - 5 для варианта 1; Рис. 12. Графики и, - 5 для варианта 2; 


однослойная схема однослойная схема 


На этих рисунках значения параметра 5 1, 2, 3, 4 соответствуют сеткам 4х4, 8х8, 16х16, 32х32. Как и в 
примере 1, линия 1 соответствует решению на базе 30110185, а линия2 — МДА. Горизонтальная линия, 


обозначенная цифрой 3, соответствует точному решению. 

Из приведенных графиков следует, что элемент, построенный по моментной схеме на сетке 32х32х2 
вариант 2, имеет относительную погрешность 4 %. 

Картины визуализации поля распределения вертикальных перемещений и, для МДА и $0185 при 


радиальной схеме разбивки '4 пластины (сетка 32х32х2) показаны соответственно на рис. 15 и 16. 
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Рис. 13. Графики и, - 5 для варианта 1; Рис. 14. Графики и, - 5 для варианта 2; 
двухслойная схема двухслойная схема 
и, м и, М 
-. 004536 
в. МИЯ _ соо 
-0,01202 — -. 003528 
== -. 003024 
-0,009016 __ 00252 
-0,006011 С | -. 002016 
С] -. 001512 
о: 7 -. 001008 
0,0 Г] -. 504Е-03 
. 410Е-06 
Рис. 15. Распределение и, (МДА) Рис. 16. Распределение и, ($0110185) 


Из приведенных рисунков видно, что несмотря на качественное совпадение картин и,, относительные 


погрешности для максимального прогиба составляют: МДА — 3,8 %; $0110185 — 71 %. Такая значительная 
погрешность при использовании 30110185 объясняется тем, что разработчики использовали для 
аппроксимации геометрии и перемещений шестиузлового КЭ «функции формы» аналогичные функциям, 
примененным для восьмиузлового элемента, т. е. без принципа «вырождения» [7]. 

Обсуждение и заключение. Построенные на базе метода двойной аппроксимации матрицы жесткостей 
объемных полилинейный конечных элементов позволяют моделировать напряженно-деформированное 
состояние строительных конструкций произвольной геометрии при различных видах внешнего воздействия. 
Принципиальным отличием предлагаемой концепции от ранее известных конечно-элементных технологий 
является то, что перемещения в данном случае задаются в глобальных координатах, а компоненты тензора 
деформаций определяются в местных, в общем случае, неортогональных осях. 

На тестовых примерах показано, что объемные конечные элементы, построенные по методу двойной 
аппроксимации, обладают устойчивой сходимостью и успешно конкурируют с элементом аналогичного типа 
50110185 вычислительного комплекса АМЗУ$ Меспашса|. 

Разработанное математическое обеспечение может быть внедрено в отечественные импортозамещающие 
программные комплексы, реализующие метод конечных элементов в форме метода перемещений. 


Список литературы 

1. дЛепюемст О.С., Тау1ог В.Г. Тре Етие ИетепЕ Мефоа, ЕШЬ еФюп. Охюга, ОК: Вийегуог-Нешетапи; 
2000. 708 р. 

2. Рама У. Нивюп. Рипаатетаб оГЕтйе Еетеп! Апайуз5. Ме\у Уок, МУ: МеОга\м НШ Сотрашез; 2004. 494 р. ОВГ: 
Ьрэ:/Амр.Киби.ас.и/ Ка]апагу/Зоигсе/ЕтИе%2Ое]етет%2О0тео4/Воок$Матенса/Рипдатета15%200%20Етще%20Е1 
степ 20 Апа[уз15., 20 Ни юоп%20(2004) раЕ (дата обращения: 15.08.2023). 

3. Рагу| Г. Гогап. А РЁ Соигуе т те Етие Ететеп! Мефоа. М№е\у Уощ, МУ: СГ Епешеенпс; 2011. 836р. 
Овг: Брз://Киби.ас.и/ОогзаРах/азег Не; /е/Месватлса/ОзааЕПе/4г паКводсв/Рагу!..ГогапАЕи$Соитзе.раЕ 
(дата обращения: 15.08.2023). 


Механика 


373 


Бйр://уезииК-доп$а.ги 


Афуапсеа Епатеетия Кезеагсй (Коз1оу-оп-Ооп). 2023;23(4):365-375. е155М 2687—1653 


4. СаПо$ А. ЕеПрра. ГигодисНоп о Етйе Еетет Мешо4а5. Вош@ег, СО: Чшуегзйу оЁ Со]огадо; 2004. 791 р. 
ОВГ: Вйрз://ушсаппаттегие!. ез.\уогаргез$.сот1/2017/01/Ёет.раЁ (дата обращения: 15.08.2023). 

5. баее4 Моауеш. Амйе Еетет! Апаруз15. Тйеогу ап Аррйсапоп ий АМ5У5. НобоКеп, №7: Ргепйсе На; 1999. 
527 р. ОВГ: Бр://Йр.детесларг.6г/зсарИпаз/ГМ7З8Ллугоз/Ешие%20Ейетеп{%20Апа|уз1з,%20ТВеогу%20ап4% 
20аррИсаНоп%20\ 1 %20АМ5У$,%20.р4Е (дата обращения: 15.08.2023). 

6. Голованов А.И., Тюленева О.Н., Шигабутдинов А.Ф. Метод конечных элементов в статике и динамике 


тонкостенных конструкций. Москва: Физматлит; 2006. 391 с. 

7. Рама У. Нийоп. Рипаатетай о{ Етие ИетепЕ Апайуя5. Ме\ Уогк, МУ: МеОтам-НШ; 2004. 505 р. ОВГ: 
Брз://\мур.Кива.ас.и/Ё”_Ка!апагу/Зоигсе/ЕтИе%20еетеп%20тео4/Воок$Митенса/Е ии атетиа15%200%20Еши 
©%20Е1етеп%20Апа[у$15,%20Ниюп%20(2004).раг (дата обращения: 15.08.2023). 

8. М1ече] Гли2 Висает, КЛалз-Лагоеп Вафе. Тре Месратсу оГ5оП45 апа Утгисигеу — Шегагстса! Модейтз апа 
ше Етийе Ейетеп! бошНоп. №е\ Уощк, МУ: Зрипеег; 2011. 597 р. 

9. ГасоБ Е1зВ, Тед ВеусНКо. А Е15Ё Соигзе т Етие ЕИетет5. НобоКеп, №]: \Пеу; 2007. 319 р. 

10. Гайджуров П.П. Конечные элементы повышенной точности для решения трехмерных задач теории 
упругости. Известия вузов. Северо-Кавказский регион. Технические науки. 2003;(1):54-57. ЧЕ: 
В рз://субе|епшка.ги/агис]е/п/копесбпуе-@етету-роуузВеппоу-юсфпозИ-Фуа-гезветуа-иебтегпуй-хаЧасв-еогй- 
иргагозИ/уме\ег (дата обращения: 15.08.2023). 

11. Чигарев А.В., Кравчук А.С., Смалюк А.Ф. АМ№Уу для инженеров. Справочное пособие. Москва: 
Машиностроение; 2004. 512 с. ОВГ: В@рз://\м\у\ гезеагспзае.пе/ргое/АКтгаусвиК/ря6ПсаНоп/262729610_АМЗУБЗ - 
Ча_ттепегоу/К$/023 1753642946 139000000/АМЗУ $-Чаттепегоу .р4Ё (дата обращения: 15.08.2023). 

12. Мадепс Е., Сиуеп Г. Тре Етйе Ейетеп! Мейо4 ап4 АррПсаноп; т Епетеетие Изуте АМ5У5. Ме\и УотК, 
МУ: брипоег; 2015. 664 р. 


ВеГегепсе$ 

1. Хлепчеулст ОС, Тау1ог ВГ. Тйе Етие Еетепт! Мейтоа, Е е@ оп. Охта, ОК: Вийегуог-Нешетапи; 
2000. 708 р. 

2. Рама У Нийоп. Еипдатета!5 ор Етие Еетет! Апа[?515. Ме\и УогКк, МУ: Тве МеСтам НШ Сотрашез; 
2004. 494 р. ОВГ: В@рз://\мр.Кпби.ас.т/Ёи_Ка|ап{агу/Зоигсе/Ете%20|етеп%2Отеоа/Воок$Матенса/Еип- 
Чатета15%200Р/20Е1Ие%20Е|етеп{%20Апа[уз1$,%20Ниюп%20(2004).раЁ (ассеззеа: 15.08.2023). 

3. Рагу! Г. Гогап. А Рй5Ё Соигуе т ше Етийе Еетет! Мефоа. Ме\м Уогк, МУ: СГ Епзшеегис; 2011. 836 р. 
ОКвГ: Брз://киби.ас.и/ОогзаРах/азеге/Ше/Месватса/ОзааЕЦе/Аг_паКВо4св/РагуП..Г.огапАЕи${Сопгзе.раЁ 
(ассеззеа: 15.08.2023). 

4. Са|оз А ЕеПрра. Гигодисйоп ю Етийе Еетет Мейо4. Во4ег, СО: Чшуегзцу оЁ Со]огадо; 2004. 791 р. 
ОВГ: В@рз://ушШсапваттегие!. ез\уогаргез$.сот/2017/01/ет.раЕ (ассеззеа: 15.08.2023). 

5. Заее4 Моаует. Рлийе Еетеп! Апа/у515. Тйеогу апа АррПсаноп ий АМ5У5. НофоКеп, №): Ргеписе НаП; 1999. 
527 р. ОВГ: Бр://Ар.Четлес арг. 6г/91зс1рНпаз/ТМ738/Глугоз/Е шие %20Еетеп{%20Апа[уз13, %20ТВеогу% 
20ап4%20аррПсайоп%20\1%20А№5У5.%20.раЕ (ассеззеа: 15.08.2023). 

6. Со]оуапоу АТ, Тушепеуа ОМ, ЗШеабитоу АЕ. Рийе Еететй Мефо4 т 5айсу апа Рупатсу ог Тит- 
ИЙ’аПеа Угисшгех. Мозсо\: Еигтай и; 2006. 391 р. (№ Клв$.) 

7. Рама У Нийоп. Рипаатетав ор Етйе Еетепт! Апа[уя15. Ме\у УотК, МУ: МеОтам-НШ; 2004. 505 р. ОВГ: 
Брз://\мур.Киеа.ас./Ёи_Ка!апагу/Зоигсе/ЕтИе%20еетет%2О0тео4/Воок$Митенса/ИЕ ип атентиа15%20020Е ши 
е%20Е1етеп%20Апа[у$15,%20Ниюп%20(2004).раЕ (ассеззеа: 15.08.2023). 

8. М1еие! Гли7 Висает, КЛаяз-Лагоеп Вае. Тре Месйатсу ор 5оПа5 апа 5гисигех — Шегагстса Мо4ейп® апа 
ше Етийе Еетеп! бошНоп. №е\ Уощк, МУ: Зрипеег; 2011. 597 р. 

9. ГасоБ Е1зВ, Тед ВеусНКо. А Е175Ё Соигзе т Етйе ЕИетет5. НобоКеп, №]: У Пеу; 2007. 319 р. 

10. Са!42Бигоу РР. Ешие Еетепз оЁ псгеазе4 Ассигасу Юг Зойуше 30 РгоМетз ог Еазисиу ТВеоту. 
Отлетзиу № ул. Мотй-Саисаяап Вевлоп. Тесйтса! 5сепсе5 Зетес. 2003;(1):54—57. ОВГ: 
В рз://субе|ешшКа ги/агае/л/Копесвпуе-@етету-роуузВеппоу-юсппозИ-Фуа-гезпетуа-мебтегиуй-тадасв-{еоги- 
иргаеоз/у\е\уег (ассеззеа: 15.08.2023). (ш Влз$.) 

11. Сшвагеу АУ, КтаусваК АЗ, Зтаучк АЕ. А№75 юг Епетеегу. Мо5сом: Мазртозиоеше; 2004. 512 р. ЧВГ: 
В@рз://уу\и\у гезеагспоае пе/ргое/АКтаусвяК/ря6НсаНоп/262729610_АМЗУЗ _Ч!а_шттепегоу/№1Кз/0Е3 1753542946139 
000000/АМЗУ $-Ча-штепегоу.раЁ (ассеззеа: 15.08.2023). (ш Виз$.) 

12. Мадепс Е, Сиуеп Г. Тре Рийе Еетепт! Мефоа ап4 АррПсапопу т Епетеетте Изуте АМУУ5. Мем УогК, 
МУ: брипоег; 2015. 664 р. 


Гайдсуров П.П. и др. Применение метода двойной аппроксимации для построения матриц жесткости объемных конечных элементов 


Поступила в редакцию 29.09.2023 
Поступила после рецензирования 31.10.2023 
Принята к публикации 18.11.2023 


Об авторах: 

Петр Павлович Гайджуров, доктор технических наук, профессор кафедры строительной механики и 
теории сооружений Донского государственного технического университета (344003, РФ, г. Ростов-на-Дону, 
пл. Гагарина, 1), ОВСТО, 5РГ\-код: 6812-9718, эрр-161@уап4ех.га 


Нина Александровна Савельева, старший преподаватель кафедры строительной механики и теории 
сооружений Донского государственного технического университета, (344003, РФ, г. Ростов-на-Дону, 
пл. Гагарина, 1), ЗРП\-код: 8437-8080, ОВСТО, птазау8 6(@лтай.га 


Заявленный вклад соавторов 

П.П. Гайджуров — постановка задачи, выбор метода решения и построение математической и 
компьютерной модели, обсуждение результатов. 

Н.А. Савельева — проведение обзора, проведение расчетов, обсуждение результатов. 


Конфликт интересов. авторы заявляют об отсутствии конфликта интересов. 


Все авторы прочитали и одобрили окончательный вариант рукописи. 


Весетуеа 29.09.2023 
Веужед 31.10.2023 
Ассерд 18.11.2023 


'АБои! ше Аифогу: 
Реег Р. Са147вигоу, Ог.5с1. (Епэ.), Ргоеззог оЁ фе Эгасвага| Месвапс$ ап4 ТВеогу оЁ Эасвгез Оерагитепт, 
Поп 5ие Тесбиса|! Ошуегзиу (1, Сагагт $4., Козюу-оп-Ооп, 344003, КЕ), ЭРПМ-соде: 6812-9718, ОВСТО, 


грр-161@уапдех.га 


Мша А. баущеуа, Зетог |есбагег оЁ фе Зиасвга| Меспатс$ ап4 ТВеогу оЁ Эгисвтез Перагитепь Роп Эа 
Тесблиса! Ошуегзиу (1, Сасайи 5а., Возюу-оп-ОРоп, 344003, КЕ), ЭРГУ-соае: 8437-8080, ОВС, 


п1пазау8 6(@тай.га 


Сеатей соттфшот5йр: 


РР Са14тВитоу: рго ет за йетепь з@есйоп о зоНоп тео4 апа БаП4те оЁ таетайса| ап сотршег п1о4е|, 
415с15510п ОР Фе гезиН$. 
МА Зауееуа: сопдисипе а геме\, сопршайопа]| апа1уз1$, 1$с15$10п оРШе гези 5. 


Соп1сЕ о 1теге5! зкиетет!: 1е алфогз 40 по Вауе апу соп1сЕ оЁ пегезй. 


АП ашфогу йауе геа4 ап4 арргоуе4 ше дпа тапизстри. 


Механика 


375 


